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Antiforromangotic (AFM) order in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) supercon- 
ducting state is analyzed on the basis of the Ginzburg-Landau theory. To examine the possible 
AFM-FFLO state in a heavy fermion superconductor CeCoIns, we focus on the incommen- 
surate AFM order characterized by the wave vector Q = Qo ± (fine with Qo = (ir, tv, it) and 
(fine || [HO] or [110] in the tetoragonal crystal structure. We formulate the two component 
Ginzburg-Landau model to discuss two denerate incommensurate AFM states with qi nc | [110] 
and [110]. Owing to the broken translation symmetry in the FFLO state, multiple phase di- 
agram of single-q phase and double-q phase is obtained under the magnetic field along [100] 
or [010]. Magnetic properties in each phase are investigated and compared with the neutron 
scattering and nuclear magnetic resonance (NMR) measurements. An ultrasonic measurement 
is proposed for a future experimental study on the possible AFM-FFLO state in CeCoIns. The 
field orientation dependence of the AFM order in CeCoIn 5 is discussed. 
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1. Introduction 

More than 40 years ago, Fulde and Ferrell, 1 ' and 
Larkin and Ovchinnikov 2 ' proposed the appearance of a 
spatially modulated phase in a spin polarized supercon- 
ductor. The original BCS theory is based on a conden- 
sate of Cooper pairs with vanishing total momentum. On 
the other hand, this Fulde-Ferrcll-Larkin-Ovchinnikov 
(FFLO) phase is composed of Cooper pairs with finite 
momentum. A spontaneous breaking of spatial symme- 
try is implied by the internal degrees of freedom of the 
FFLO phase arising from inversion or reflection symme- 
try. 

The search for the FFLO phase has been pursued since 
its proposal with mixed success. Indeed it has been dis- 
cussed for various systems such as superconductors, 3 ~ 9 ' 
cold atom gases, 10 ~ 12 ' and quark matter. 13 ) Nevertheless, 
the discovery of a new superconducting phase in CeCoIn 5 
at high magnetic fields and low temperatures 3,4 ' came as 
a surprise and triggered much interest, as it satisfied all 
the immediate criteria for an FFLO state, sitting in the 
right parameter space in the H-T-pha.se diagram. 14 ' 

The high-field superconducting (HFSC) phase of 
CeCoIn 5 has been examined based on the FFLO- 
concepts by many groups. 14 ~ 23 ' Recent observations of 
magnetic order coexisting with the HFSC phase point 
towards a more complex situation and demands the reex- 
amination and extension of the basic picture. 24 ' 25 -' Neu- 
tron scattering measurements show that the wave vec- 
tor of AFM order is incommensurate, Q = Qo ± <fi nc 
with Qo = (7r,7r,7r) and q inc ~ (0.12tt, ±0.12tt, 0), 
and independent of the orientation of inplane magnetic 
field. 25,26 ' The ordered Ce-moments Maf are oriented 
along the c-axis. For the magnetic field along [110], Qinc 
is perpendicular to the applied magnetic field <fi nc ~ 
(0.12-7T, — 0.12-7T, 0). The data for the magnetic field along 
[100] direction may be best interpreted as being due to 



two degenerate incommensurate AFM states with $ nc ~ 
(0.127r,0.127r,0) and q inc ~ (0.12tt, -0.12tt, 0), 25,26 ' as 
we will discuss below. 

We may conclude that this magnetic order is stabilized 
by superconductivity because it is absent in the normal 
state. This strong cooperation between the magnetism 
and superconductivity contrasts the behavior of other 
heavy fermion superconductors where the magnetic order 
competes usually with superconductivity. 27 ' It is worth 
noting that the HFSC phase is enlarged by pressure. 19 ' 
This feature is also different from the other magnetic 
phase diagram in Ce-based heavy fermion systems. 27 ' 

The observation of this unconventional magnetic order 
initiated a number of attempts for theoretical explana- 
tions. 28 ~ 36 ' Here we propose a theory based on the pres- 
ence of an AFM quantum critical point near the super- 
conducting phase of CeCoIn5. 37 ~ 41 ' As we have shown 
previously, AFM order appears within the inhomogc- 
neous Larkin-Ovchinnikov state in the vicinity of the 
AFM quantum critical point. 28,29 ' The mechanism cou- 
pling magnetism and FFLO superconductivity is consis- 
tent with the fact that the magnetic order is restricted 
to the superconducting phase in the H-T phase dia- 
gram. 24,25 ' Alternative proposals are based on the emer- 
gence of a pair density wave state together with magnetic 
order in the HFSC phase. 34 " 36 ' 

In order to identify the HFSC phase of CeCoIn 5 it 
is important to find unambiguous ways to compare the 
proposed phase with experimental findings. Here we in- 
tend to examine the AFM-FFLO state proposed by us 
previously. First, we show that several phases can ap- 
pear in the AFM-FFLO state when the magnetic field is 
applied along [100] or [010] direction. Second, we clarify 
the magnetic structure of each phases and discuss the 
consistency with the experimental results of CeCoIn 5 . 
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2. Ginzburg-Landau theory 

Our theoretical analysis is based on a phenomenolog- 
ical Ginzburg-Landau model formulated on the basis of 
the microscopic calculation in Refs. 29 and 42. and moti- 
vated by recent neutron scattering measurements of Ken- 
zelmann et aL 25 ' 26 ) We describe the AFM order in the in- 
homogeneous Larkin-Ovchinnikov state by means of the 
Ginzburg-Landau functional of the free energy 

+c 1 H x H y (r ] 2 - ri) 

~\mV2^2 c 2 {n)5{q 1 - q 2 ,2nq FFLO ). (1) 

n 

We use an order parameter with two components, r\\ and 
7/2 corresponding to the two degenerate AFM states with 

Q <~ Qo±9i„c^ an d Q ~ Qo±q^nJ ■ The wave vectors of in- 
commensurate AFM state in the system with full transla- 
tional symmetry are given by q^J = (0.1257T, 0.1257T, 0) 
and q^J = (-0.125tt, 0.125tt, 0). In our study the in- 
commensurate wave vectors q^J and q^J arc not mi- 
croscopically derived but assumed on the basis of the 
experimental results in Refs. 25 and 26. We think that 
the q^J and q^J are pinned through nesting features 
in the band structures. Note that the order parameter 
(771,772) is not a vector, but a director. 

We describe the magnetic moment M(r) = 
e % Q°' r M&F{r) where M\ F {f) is a slowly varying ampli- 
tude of the magnetic moment perpendicular to the ap- 
plied magnetic field at f = (x, y, z). Then, the AFM stag- 
gered moment is given by 

M AF (r) = M [r]i cos(<fi • r ) + t/ 2 cos(<f 2 ■ r)]. (2) 

Both M and F are scaling factors to write the 
Ginzburg-Landau free energy in a dimcnsionless form 
(eq.(l)). Moreover, is the "bare" Neel temperature 
(for ci = c 2 (n) = 0, i.e. unrenormalized by coupling to 
the magnetic field and to the FFLO-modulation of the 
superconducting phase) , and £af is the correlation length 
of the AFM order. The inplane lattice constant is chosen 
as the unit length. 

Choosing the coupling constant b positive we stabi- 
lize the "single-g" magnetic structure (7/1,772) oc (1,0) 
or (0,1), while a negative value for b would favor the 
"double-<7" magnetic structure (771,772) oc (1,1). From the 
microscopic point of view 6 > is more appropriate, since 
the single-q magnetic phase gains more condensation en- 
ergy than the double-g phase in most cases. We assume 
b > in this paper, while the cases of negative b have 
been discussed in Ref. 43. 

The constant c\ denotes the coupling strength of the 
order parameter to the external magnetic field H = 
(H x ,H yi H z ). Comparing with the neutron scattering 
data for H | [110], 25 ' c\ must be positive in CeCoIn 5 , 



as to pick the incommensurate wave vector q^J _L [110]. 
This choice is consistent with our theoretical analysis for 
the AFM-FFLO state, 28 - 29 ' 42) while it is incompatible 
with the theoretical scenario without taking the FFLO 
state into account. 32 ' 33 ) When we consider the magnetic 
field along the [100] (or [010]) direction, this term is inac- 
tive and, therefore, the degeneracy of 7/1 and 772 remains. 
We focus on this situation in this paper. 

The effects of the broken translational symmetry in the 
inhomogeneous Larkin-Ovchinnikov state are taken into 
account in the last term (commensurate term) of eq.(l). 
We define 5{q\ — q 2 , 2ti<7fflo) = 1 when the commensu- 
rate condition q\ — q 2 = 2nq FF Lo is satisfied for an integer 
77 and otherwise 5{q\ — q 2 , 27m/fflo) = 0. The modula- 
tion vector of FFLO state is denoted as <7fflo , and then 
the order parameter of superconductivity is described 
as A(f) = A sin((fFFLO • r)- The commensurate term 
describes the pinning effect of FFLO nodal planes for 
the AFM moment. According to the microscopic calcu- 
lation based on the Bogoliubov-de-Gennes (BdG) equa- 
tion, the AFM moment is enhanced around the FFLO 
nodal planes where the superconducting order parameter 
vanishes. 28 ' 29 ' 42 ^ Then, we obtain the positive coupling 
constant c 2 (n) > 0. 

3. Phase Diagram 

Let us now consider the FFLO state in a case with 
2ng*FFLO is close to q^J — q^J for one integer n = N, 
i.e. we assume q^J — q^J <~ 2A<fpFLO- This has conse- 
quences for the structure of the magnetic order. When 
the wave vector q^J — q^J is not in the vicinity of 
27i<7fflo f° r an Y integer n, the single-q magnetic phase 
with (771,7/2) oc (1,0) or (0,1) is realized. 

We determine the magnetic structure by minimizing 
the Ginzburg-Landau free energy (eq.(l)) with respect 
to the order parameters 7/1 and rj 2 and their momen- 
tum qi and q 2 for given temperature T and FFLO wave 
vector (/fflo • We assume 771 > r/ 2 > without loss of 
the generality. As for the direction of qpFLO , we assume 
9fflo = <7fflo^ for H | [100] as in Refs. 29 and 42. 
Then, the phase diagram is determined by the renor- 
malized parameters T/T^ and £af<7o, where go is de- 
fined by q^J - q^J + 2q x = 2Nq FFho . The parameter 
qo describes the mismatch of the harmonic FFLO wave 
vector 27V<2fflo and the incommensurability along the 
a;-axis q^J — q^J ■ Note that q a = when the condi- 
tion q-^2 — q^J — 2Nq FFLO is satisfied. Since the FFLO 
wave number (/fflo grows with increasing magnetic field 
and/or decreasing temperature, 14 ) qo increases with the 
magnetic field, and changes its sign at the line in the 
H-T phase diagram. 

We find two possible phases. One is the single-g phase 
where 

771 > 7/2, (3) 

Qi = QinJ + i 1 - x )lox, (4) 
92 = QinJ - (1 + x)q x. (5) 
Note that x lies between and 1 and is determined by 
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minimizing the free energy. The other is the double-q 
phase where 



m = V2, 



Qi 



<12 



= q W 

^inc 
,(2) 



q x. 



(6) 
(7) 
(8) 



By analyzing the quadratic terms of eq.(l), we find that 
the double-g phase is stabilized immediately below the 
Neel temperature for (£af<?o) 2 < C2(iV)/8, while a single- 
q phase is stabilized for (£af<7o) 2 > C2(N)/8. 
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Fig. 1. (Color online) (a) Phase diagram of Ginzburg-Landau 
model in eq.(l) for £af<70 and the renormalized temperature 
T/T^J. The definition of qo is given in the text. The go increases 
with the magnetic field. The inset shows the same phase diagram 
scaled up around the Neel temperature, (b) The ratio of order 
parameters 772 jr\\. 



Ref. 43. We here focus on the case of Fig. 1 and examine 
the consistency with experiments. 

Figure 1 (b) shows the ratio of order parameters r\\ and 
r]2- We find that 772/771 — 1 m the double-q phase while 
the ratio decreases in the single-q phase with decreas- 
ing temperature T/T^ and/or increasing the mismatch 
I£af<7o|- 
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Fig. 2. (Color online) The £af<?0 dependences of order param- 
eters r)i and r]2 at (a) T = 0.5T°, (b) T = 0.9T°, and (c) 
T = 0.99T£}, respectively, (d) Temperature dependence of r]i 
and r]2 at £af<30 = 0. 



The behavior of the order parameter can be seen in 
the plot of the £af<Zo dependence of order parameters for 
T = 0.5T£, T = 0.9T$, and T = 0.991$ in Figs. 2(a)- 
2(c), respectively while the temperature dependence for 
£af<?o = is shown in Fig. 2(d). The phase transition 
between the doublc-q phase and single-q phase is a con- 
tinuous second-order transition. This phase transition is 
characterized by the broken mirror symmetry with re- 
spect to the x- and y-axes in the single-q phase. 



Figure 1(a) shows the phase diagram obtained by the 
numerical minimization of the free energy in eq.(l) using 
the parameters b = 0.1 and c 2 (N) = 0.01. The single-q 
phase is stable in most parameter regimes, because the 
quartic term br^rfa dominates. Only near the Neel tem- 
perature quadratic terms are leading and can stabilize 
the double-q phase for small |£af9o| < \/c2(N)/8. As 
shown in the inset of Fig. 1(a), the Neel temperature in- 
creases slightly around the commensurate line £af<7o = 
due to the locking-in effect of the magnetic and the FFLO 
modulation. However, the enhancement of the Neel tem- 
perature Tn — = 62 ^ Tjj is rather small, unless the 
coupling constant C2(N) is large. Various phase diagrams 
for several parameter sets (b,c(N)) have been shown in 




XXX 



Fig. 3. (Color online) The spatial dependences of AFM staggered 
moment M\p(r) at r = (x,y, 0). (a) Single-q phase at T = 
0.5T°, (b) Single-? phase at T = 0.9T°, and (c) Double-g phase 
at T = 0.95T^, respectively. We choose £af<?0 = 0. 



We show the spatial dependences of the AFM stag- 
gered moment Maf(t^ for several parameters in Fig. 3. 
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Since the AFM moment is uniform along the z-direction, 
we show the Mxeif) in the x-y plane. The checker- 
board magnetic structure is realized in the double-^ 
phase (Fig. 3(c)), while it changes to the stripe magnetic 
structure in the single-g phase (Figs. 3(a) and 3(b)). We 
investigate these phases in more details and discuss the 
experimental results of CeCoIns in the next section. 

4. HFSC Phase in CeCoIn 5 

In order to discuss the possible AFM-FFLO states in 
CeCoIns on the basis of our Ginzburg-Landau analysis, 
we show the schematic phase diagram of CeCoIns for the 
magnetic field and temperature in Fig. 4. The amplitude 
of FFLO modulation vector qfflo has a maximal pos- 
sible value of (?fflo ~ V£ with £ being the coherence 
length of the superconducting state. According to the 
experimental estimate of £, the minimum number of N, 
which satisfies the condition q^J — q^J = 2^Vg*FFLO for 
commensurate wave vectors, is approximately N ~ 4. 
This condition is satisfied in the FFLO state on a se- 
quence of commensurate lines (dashed lines in Fig. 4) 
with N = 4,5,6,7...., and the double-q phase is stabi- 
lized around the AFM transition line (shaded area of 
Fig. 4). Note that the double-^ phase does not appear 
close to the first-order normal-to-FFLO transition line at 
which the AFM moment as well as the superconducting 
order parameter appear discontinuously. Since the cou- 
pling constant C2(N) decreases with increasing N, the 
double-g phase is suppressed in the low magnetic field 
region. 



AFM-FFLO 




surate lines can appear in the phase diagram and stabi- 
lize the double-g phase around T = Tj$. However, it is 
expected that the effect of the vortex lattice on the mag- 
netic order is smaller than that of the FFLO nodal planes 
when the Maki parameter is large enough to make the 
lattice spacing of vortices much larger than the coherence 
length. 

4-1 Neutron scattering 

The neutron scattering measurements have well de- 
termined the structure of AFM order and its magnetic 
field dependence. 25 ' 26 ^ For the magnetic field along [100] 
direction the elastic Bragg peaks appear at both Q = 
Qo±q i i 1 J and the symmetry related Q = Q ±q i £ ) - 25,26) 
At first sight, this result seems to be incompatible with 
the single-g phase in which, for instance, the Bragg peaks 
at Q = Qo ± q inc should be much weaker than those 
at Q — Q ± q^J (see Fig. 5(a)). However, taking into 
account that the sets of wave vectors correspond to de- 
generate single-g phases in a field along [100], it is rather 
likely that the domains have formed of the two states 
|?7i| > l r /2| and \r]i\ < 1 772 1 - Because the domain struc- 
ture is history-dependent, it can be identified in future 
experiments, for example, by changing the field orienta- 
tion slightly away from [100] direction in order to lift the 
degeneracy. 
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Fig. 5. (Color online) The weight of Bragg peaks defined by 
\M(q)\ 2 with M(g) = J2 P M{r)e^ and q = ( 9 x,<Zy,7r). The 
summation is taken over the Nq lattice sites, (a) Single- 
q phase at T = 0.5T°, (b) Single-g phase at T = 0.9T°, and 
(c) Double-g phase at T = 0.95T^J, respectively. We choose 
?AF90 = 0. 



T 

Fig. 4. (Color online) Schematic phase diagram for the magnetic 
field along [100] direction. "BCS", "FFLO", and "AFM-FFLO" 
states are shown in the figure. The dashed lines show the com- 
mensurate lines which are explained in the text. The shaded 
region shows the double-g phase in the AFM-FFLO state. 



It should be noticed that the AFM-FFLO phase is 
mostly covered by the single-g phase. Therefore, the ex- 
perimental results of CeCoIns should be compared with 
the properties of the single-g phase. There are ways to 
distinguish the single-q phase experimentally from the 
double-g phase. Indeed the experimental results are con- 
sistent with the single-g phase, as we will discuss below. 

When we take into account the broken translational 
symmetry arising from the vortex lattice, other commen- 



As discussed above and stated in eqs.(4), (5), (7) and 
(8), the incommensurate wave vectors q± and q\ depend 
on the magnetic field. However, the deviation from q^J 
is rather small for the main Bragg peak, with \q\— q^J \ < 
£aW c 2( N )/ 8 - For the parameters c 2 (N) = 0.01 and 
£af = 3, we obtain \q\ — q^ c '\ < 0.0047T. Moreover, the 
shift of <fi rapidly decreases upon lowering the tempera- 
ture below Tn (see Fig. 6). Therefore, additionally chal- 
lenged by the tininess of the magnetic moment near the 
Neel temperature, most likely the shift of Bragg peaks is 
experimentally unobservable. Indeed, in the experiments 
no magnetic field dependence of <fi has been observed so 
f ar 25,26) Nevertheless, it would be interesting to search 
for a possible shift of q\ in future experiments. 

Finally, we would like to state that also the case of 
b < stabilizing the double-g phase shows four Bragg 
peaks consistent with the neutron scattering data for 
H || [100] without having to assume domain formation. 
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Fig. 6. (Color online) The shift of the position of main Bragg 
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peak from qS^ . We plot £j±pAq for T ■■ 



0.5Tg 



: 0.99T°, 



where Aq is defined by Aqx = qi — q. 



with two peaks at the upper and lower edge of the dis- 
tribution, if (771,7/2) <x (1,0) or (0,1) (single-q phase as 
in Fig. 7(a)). A single peak in the center of the distribu- 
tion (h = 0) is found for \rjxl = |t7 2 | (double-g phase as 
in Fig. 7(c)). In Fig. 7 we show the evolution of the dis- 
tribution function P(h) for the same parameters as used 
for Figs. 3(a-c). 

The experimentally observed Knight shift distribution 
has the shape of a broad double peak structure. This 
is consistent with the realization of a single-q phase and 
clearly does not fit to the double-q phase. Thus, the phase 
diagram with the AFM-FFLO state, which is almost en- 
tirely covered by the single-g phase, is consistent with 
the actual NMR measurements. 



However, this case yields a sizable shift of q\ at low tem- 
perature, and therefore seems to be incompatible with 
the neutron scattering data. 43 - 1 Moreover, as we will elu- 
cidate now this scenario is incompatible with the NMR 
measurements. 

4.2 NMR 

The first evidence for AFM order in the HFSC phase 
of CeCoIns has been obtained by the NMR measure- 
ment. 24 > The analysis of the NMR spectrum pointed to- 
wards incommensurate AFM order. This interpretation 
has subsequently been confirmed by the neutron scatter- 
ing measurements in a decisive way. 25,26 ) Recent analysis 
showed that the direction of AFM moment along c-axis is 
also consistent with the NMR data by assuming dipolar 
hyperfine coupling. 44 ~ 47 ) 

Applying the magnetic field parallel to the a;-axis there 
are three distinct In-sites for NMR, the inplane In(l)-site 
and the out-of-plane In(2a)- and In(2b)-site. The latter 
two lie on yz- and zx-planes, respectively, relative to 
the Ce-layers. Looking at the dipolar fields due to the 
magnetic moments on the Ce-sites only the In(2b)-site 
develops a field component parallel to the x-axis which 
contributes to the Knight shift. Restricting to the two 
nearest Ce-ions of a In(2b)-site at r we obtain for the 
dipolar field, 

H x {r) = A[M(r + b+) - M(r + &_)], (9) 

with A a constant and b± = (±1/2, 0, d) with d the dis- 
tance of In(2b) site from the Ce-plane. Note that there 
are In(2b)-sites above and below the Ce-plane. Using 
eq.(2) we derive, 

H x (f) = 

A'M cos(Q • r) [m cos(qi • r) + n 2 cos(g 2 ■ r)} , (10) 

which eventually leads to the field distribution P(h) ob- 
servable in the Knight shift, 

r 

The sum runs over Nq sites of square lattice representing 
the In(2b)-sites. 

It is obvious that P(h) shows a broad distributions 



(a) T=0.5T N ° (b) T=0.9T N ° (c) T=0.95T, 




Fig. 7. The distribution function of the internal field P(h) at the 
In(2b) sites given by eq.(ll). The parameters are the same as in 
Figs. 3(a), 3(b) and 3(c), respectively. We choose AMq = |. 



We now turn to the relationship between the analysis 
given in this paper and our previous microscopic study 
of the spatial structure of AFM-FFLO state. We have 
shown that three basic forms of the AFM-FFLO state 
are possible. 42 ) We distinguish (1) the 'extended' case 
for which the magnetic moments are slightly suppressed 
around the nodal planes of FFLO-superconducting or- 
der parameter, (2) the 'weakly localized' case in which 
the magnetic moments are slightly enhanced around the 
nodal planes, and (3) the 'strongly localized' case with 
a magnetic order parameter very localized around the 
nodal planes. Our phenomenological theory in this pa- 
per is justified in the cases (1) and (2), while it is invalid 
for the case (3) because the incommensurate wave vec- 
tor <fr nc would not be parallel to [110] and [110] for the 
magnetic field along [100] axis. According to our analysis, 
only the case (2) is consistent with the neutron scattering 
measurements in which the position of magnetic Bragg 
peaks is independent of the magnetic field and g! nc is 
perpendicular to the field (q\ nc _L H) for H \\ [110] or 
H || [110]. 25 ' 26 ) 

In the case (2), a weak enhancement of magnetic mo- 
ment around the FFLO nodal planes, which is not taken 
into account in this paper, affects the NMR spectrum. 
This results in tails at the edges of the distribution func- 
tion P(h). This additional contribution to P(h) does, 
however, not alter the qualitative feature of the double 
peak structure anticipated for the single-g phase. 

We here discuss the NMR measurements at In(l)- 
and In(2a)-sites. Since the AFM moment does not (or 
only weakly) yield internal field at these sites, we may 
avoid the effect of AFM order and elucidate the magnetic 
properties of FFLO state without AFM order. Recent 
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NMR measurement at In(l)- and In(2a)-sites actually 
reports an evidence for the FFLO superconducting state 
in CeCoIn 5 47 ) 

4-3 Ultrasonic measurement 

As we have shown in Fig. 4, the region of the AFM- 
FFLO state of CeCoIn5 is dominantly the single-g phase. 
Because the order of single-q phase has Zi character, a 
probe sensitive to the crystal symmetry could give fur- 
ther helpful evidence for this phase. Ultrasound measure- 
ments provide such an experimental tool, as we discuss 
here. 

AFM order and crystal lattice deformations are cou- 
pled to each other through the band structure and spin- 
orbit coupling. Avoiding microscopic details, we may de- 
scribe this interplay by means of a general Ginzburg- 
Landau free energy, 

F> = F( m ,r) 2 )+F el _ ai + F eh (12) 

^el-af = iTiOxx + eyy) + ^zzjivl + vl) 

+73£ xy (?7i -ill), ( 13 ) 

1, 



-[Cn(e xx + £y y 



+ C33£zz 
CxxCyy + 2Cl3(e xx + £yy)£zz 

+4C 44 (e yz + e xz )+4C 66 e xy ], 



yz 1 xz/ 1 Du '-xyJ 3 

(14) 

where the coefficients represent elastic constants of 
a tetragonal crystal lattice, and = \{duijdxj + 
duj/dxi) defines the strain tensor with u being the local 
displacement vector. 7, are coupling constants between 
strain and magnetic order parameter. 

We examine now the coupling to ultrasound modes, 
longitudinal (L) and transversal (T1,T2) with inplanc 
and c-axis propagation directions. For the propagation 
along [100] the following strain fields are involved, 



L: e. 



Tl: e 



xyi 



while for [110], 



L. €xx ~~ £yy ~~ f-xyi Tl. E x 



T2: c xzi 



^yyi T2. € xz 



and for [001], 



T. € XZ , €y 



(15) 



(16) 



(17) 



Note that for Tl (T2) the transverse polarization lies in 
the basal plane (along the z-axis) and for T it has two 
independent polarization directions in the basal plane. 

The sound velocity corresponding to these modes can 
be renormalized at the phase transition to the AFM- 
FFLO phase. This corresponds to the renormalization of 
the elastic constants, given by 

d 2 F' „ ^ 8 2 F' drj n 



Co 



dt n der- 



71=1,2 



de a dn n det 



(18) 



at equilibrium, i.e. 

dF' _ dF' 
dr] n de a 

We find coupling to the ultrasound mode, if the corre- 



. 



(19) 



Table I. Coupling of ultrasound modes to single- and doublc-g 
phases: "yes" or "no" denote which elastic constants are renor- 
malized or unrcnormalizcd. 



mode 


single- q 


double-g 




r-i nnl 








T 

l_j 


yes 


yes 


I'll 


Tl 


yes 


yes 


Cm 


T2 


no 


no 


C44 


[110] 








L 


yes 


yes 


C11 + C12 + 2Cee 


Tl 


no 


no 


C11 ~ C12 


T2 


no 


no 


C44 


[001] 








L 


yes 


yes 


C33 


T 


no 


no 


C44 



Table II. Anomaly of ultrasonic sound mode at the magnetic 
transition. The sound velocity (and elastic constant) shows a 
"kink" , "jump" , or "divergence" . In case of the double transi- 
tion, magnetic transition to the double-q state occurs at T = 
and that to the single-q state occurs at T = T^2- In case of the 
single transition, the transition to the singlc-g state occurs at 
T = T N . 





double transition 


single transition 


mode 


T = T N T = T N2 


T = T N 


Tl for [100] 


kink divergence 


jump 


L for all 


jump jump 


jump 



sponding terms, 



d^el-af 
de a dn n 



(20) 

«q=0;jj„=jj„o 

are finite. We give now a list of couplings for the two 
phases, single-g [(771, 772) oc (1, 0)] and double-q [(771, 772) oc 

(1,1)]- 

Table I shows that both single- and double-q phases 
couple in the same pattern to the different modes. Thus 
there are no selection rules distinguishing the two phases. 
Important is, however, to notice that both phases couple 
not only to longitudinal modes, but also to a transverse 
(Tl) mode propagating in the [100] direction. This is a 
signature of the multi-component (magnetic) order pa- 
rameter, since a single component magnetic order, such 
as the commensurate AFM order with Q = Q , does 
not couple to any transverse mode. The anomaly in the 
Tl mode arises from the last term of eq.(13) which is 
described as 73£ xy Q xy with use of the quadrupole order 
parameter Q xy — n\ — n\. Thus, the internal degree of 
freedom having a quadrupole symmetry manifests itself 
in the transverse sound mode, as in the non-magnetic 
quadrupole order appearing in many heavy fermion sys- 
tems. 48 ) 

We find characteristic behaviors of this sound mode 
and summarize in Table II. For |£af9oI < ^/c 2 (A)/8 
the double magnetic transitions occur at T = Tn and 
T = Tn2 (Tn2 < Tn). Then, the Tl sound mode prop- 
agating along [100] direction shows a kink at T = Tn 
and diverges at T = Tn2- On the other hand, when the 
single magnetic transition to the single-q phase occurs 
for I^af^oI > ^/c 2 (A)/8, this sound mode shows a jump 
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at T = Tn- Thus, these two cases can be distinguished 
by the ultrasonic measurement for the Tl mode along 
[100] direction. All longitudinal modes show a jump at 
the magnetic phase transition. 

The ultrasonic measurement by Watanabe et al. actu- 
ally found the anomaly in the Tl mode along [100] di- 
rection, but it has been attributed to the modulation of 
superconducting order parameter in the FFLO state. 15 ) 
More detailed ultrasonic measurement with focus on the 
AFM order is highly desired to identify the HFSC phase 
of CeCoIn 5 . 

It is obvious from eq.(13) that a uniaxial stress along 
[110] or [110] direction (e xy ^ 0) would split the degen- 
eracy between iji and 772 such that a single domain phase 
could be obtained for the single-q phase. 

5. Summary and Discussion 

In this study we investigated the coexistence of the in- 
commensurate AFM order with the FFLO superconduct- 
ing state to interpret the properties of the HFSC phase 
of CeCoIn 5 . We assume that two incommensurate AFM 
modulations with <f; nc || [110] and <7; nc || [110] are degen- 
erate for the magnetic field along [100] direction. Taking 
the coupling of the magnetic phase and the FFLO mod- 
ulation of the superconducting phase into account we 
show that there could be multiple AFM phases in the 
H-T phase diagram. The so-called "single-g" phase cor- 
responding to one of the two incommensurate wave vec- 
tors (fi nc is dominant in the phase diagram. However, the 
"double-<?" phase superposing both wave vectors can be 
induced under certain conditions, if there is commensu- 
ration condition between FFLO and AFM modulations. 
This phase is unstable against the second order transition 
towards a single-g state as our Ginzburg-Landau model 
demonstrates. 

We have shown that the single-g phase is consistent 
with the experimental results of the NMR and neutron 
scattering measurements by taking into account the do- 
main structure of two degenerate single-g phases. A fu- 
ture experiment by the ultrasonic measurement has been 
discussed. 

Finally, we discuss the field orientation dependence of 
the AFM order in CeCoIns . Since the staggered magnetic 
moment lies along the c-axis owing to the spin-orbit cou- 
pling, the AFM order should be suppressed by the mag- 
netic field along [001] direction. Furthermore, a small H C 2 
along c-axis indicates a weak paramagnetic effect which 
is unfavorable for the magnetic order. 31 ' 32 ) Thus, it is ex- 
pected that the AFM order does not occur for this field 
direction even when the FFLO superconducting state is 
stabilized. Indeed, recent experiments have shown that 
the magnetic order is suppressed for H || [001] , 49 ' 50 ) 
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